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I present a theoretical treatment of parametric scattering 
in strong coupling semiconductor microcavities to model ex- 
periments in which parametric oscillator behaviour has been 
observed. The model consists of a non-linear excitonic os- 
cillator coupled to a cavity mode which is driven by the ex- 
ternal fields, and predicts the output power, below threshold 
gain and spectral blue shifts of the parametric oscillator. The 
predictions are found to be in excellent agreement with the 
experimental data 

Recent experimental studies have demonstrated that 
very large optical non-linearities can be obtained in res- 
onantly pumped strong-coupling microcavities. The ex- 
citations of these structures are polaritons, mixed modes 
which arc part exciton, part cavity photon, and the non- 
linearity is due to interactions between the exciton com- 
ponents, which cause the polaritons to scatter off each 
other. This leads to a parametric process where a pair of 
pump polaritons scatter into non-degenerate signal and 
idler modes, while conserving energy and momentum. 
The scattering is particularly strong in microcavities be- 
cause the unusual shape of the dispersion, shown in the 
inset to Figjl], makes it possible for pump, signal and 
idler all to be on resonance at the same time. 

A further important property of planar microcavities 
is the correspondence between the in-plane momentum 
of each polariton mode and the direction of the external 
photon to which it couples. This makes it quite straight 
forward to investigate parametric scattering using mea- 
surements at different angles to access the various modes. 
Two types of process have been studied in this way: para- 
metric amplification, where the scattering is stimulated 
by excitation of the signal mode with a weak probe field, 
and parametric oscillation, where there is no probe and 
a coherent population in the signal mode appears spon- 
taneously. 

Parametric amplification in microcavities was first ob- 
served by Savvidis et alu, using ultrafast pump-probe 
measurements. The structure was pumped on the lower 
polariton branch at an incident angle of 16.5°. Narrow 
band gains of up to 70 were observed in the region of 
the polariton feature for a probe at 0°, along with idler 
emission at 35°. These arc a set of angles for which the 
pair scattering resonance condition is satisfied. A related 
scattering process has been studied by Huang et afl us- 
ing two pump beams at ±45° and a probe at normal 
incidence. The experimental results of Ref. [f] have been 
modelled by Ciuti et am using a microscopic quantum 
treatment of polariton-polariton interactions. 

Parametric oscillator behaviour has been observed 
Stevenson et ala and Baumberg et alB, in CW exper- 



iments with the pump incident on the lower polariton 
branch at the 'magic' angle of about 16°. Above a thresh- 
old pump intensity, strong signal and idler beams were 
observed at about 0° and 35°, without any probe stimu- 
lation. The coherence of these beams was demonstrated 
by significant spectral narrowing, proving that they are 
due to a parametric process rather than resonantly eft- 
hanced incoherent photoluminescence. Houdre et aM, 
have also observed a nonlinear emission at 0° for a struc- 
ture pumped at 10°. However, in this experiment the pair 
scattering resonance condition is not satisfied, suggesting 
that different physics may be involved. 

The purpose of this paper is to develop a simple clas- 
sical model which provides a unified treatment of both 
amplifier and oscillator in the CW regime. This model 
is based on the textbook treatment _pf parametric phe- 
nomena in systems such as LiNbOaO. Indeed, the mi- 
crocavity behaviour has similar characteristics to a typi- 
cal doubly-resonant parametric oscillator, where just the 
signal and idler modes are cavity resonances - the pump 
resonance is mainly important in enhancing the strength 
of the nonlinear effects. However, there are also signif- 
icant novel aspects to the model: the microcavity op- 
erates in the strong coupling regime, where the modes 
are cavity polaritons not simple photons, and instead of 
a non-resonant x nonlincarity, the exciton provides a 
highly resonant x effect. 

I. MODEL 

The theoretical model is a classical treatment of a non- 
dispersive exciton mode, ip(r), with energy lo x , coupled 
to a cavity mode, </>(r), with dispersion w c (fc), which is 
driven by the external fields. To account for broaden- 
ing processes, ui x and u> c (k) are taken to be complex 
energies, with imaginary parts "f x and 7 C respectively. 
The exciton mode is non-linear, with potential energy 

To model the parametric processes, the cavity is driven 
by harmonic plane waves, consisting of a pump with am- 
plitude F p at (o! p ,kp) and a probe with amplitude F s at 
(w s ,k s ). The cavity and exciton modes are also expressed 
as a sum of plane waves at (w p ,k p ) and (w s ,k s ), plus an 
idler at (w, = 2u p — lu s , k; = 2k p — k s ). The cavity mode 
is linear, so the equations of motion for the pump, signal 
and idler mode separate out, giving 

(uJc(kp) 2 - u p 2 )<j)p + gtpp = Fp (1) 
(uJ c (k s ) 2 - LO s 2 )(f) s + gip s = F s 
(w c (fc l ) 2 - uji 2 )4>i + gi>i = 0, 



1 



where g is the strength of the coupling between the exci- 
ton and the cavity photon. The exciton equations are 
more complicated because the non-linearity generates 
many terms at different frequencies and wave-vectors. 



Only the terms at frequencies u>. 



are retained 



effectively ignored, though it can be considered to be in- 
cluded as a renormalisation of the polariton frequencies: 
at each point on the dispersion, it leads to a blue shift of 
approximately 2k\ x\ \ip P \ 2 - For low pump powers, where 
the blue shift is small compared to the pump polariton 



here: the others are at very different frequencies, such as width, and when ui p 



3ui p , or are weak, less than 0(ip p ). This leaves 



ing Eq.(ky), by tp p 



ujp, ipp can be approximated, us- 
—iCpXp* fp/"f p , and the blue shift 



{ui x 2 - ui p 2 )ip p + g<j) p + K\ipp\ 2 ip p + 2nip s ipiip p * = 
{ui x 2 - ui s 2 )tp s + g<p s + 2n\ip p \ 2 ip s + ntpp 2 ^* = 
(ui x 2 - ui 2 )ipi + g<f>i + 2n\ipp\ 2 ipi + Kip p 2 ip s * = 



(2) 



5uj° ss 2Jt|jr| 



2 \Cp\ \X p 
IP 2 



(5) 



These equations can be simplified by using Eqs.([y) to 
eliminate the cavity photon fields (p and write everything 
in terms of the exciton fields ip. It is also convenient to 



where I p = \f p \ 2 is the pump intensity. 



approximate u> x ' 
fl = g/ui x , 7 
become 



2u> x {w x — ujp) etc and define 



\k/w x and / = \F/l 



UJ X + K.\lp p \ — UJ p 



(n/2f 



ui c {kp) 



UJ X + 2k|"0 p | — LU S — 



uj x + 2k|V>p| —oJi — 



(fi/2) 2 



Ipp + 2 K1p s 1pi1p p * — 



Ips + K1pp 2 1pi* = 



Then Eqs 

-0/2 
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It is interesting to compare the present classical model 
with the treatment in Ref. @, which gives a good fit to 
the pump-probe parametric amplifier experiments of Ref. 
ED. This treatment was based on a quantum mechanical 
picture of the exciton-exciton scattering process. How- 
ever, with the approximations that were made, Eqs(l-3) 
of Ref. ^, contain essentially the same physics as Eqs. (Q) 
f p (3a) here. Of course, using a microscopic model gives a value 
for the non-linearity k. However, k only imposes a scale 
on the problem: rescaling all the fields so ip — > ip/yffi, 
f — > f/y/K, effectively makes 75=1. 



fs (3b) 



UJc(ki) 



ipi + ntpp 2 ip a * = (3c) 



Eqs.([3|) constitute the basic model for parametric pro- 
cesses in a microcavity. The terms in \ip p \ 2 represent the 
renormalisation of the exciton energy due to the pump 
population. The other non-linear terms provide the scat- 
tering, which is the main interest here: Kip p 2 ipi* and 
Kip p 2 ip s * in ( |§b| , |3c| ) describe the build up of the popu- 
lation in the signal and idler modes, while 2~Kip s ipiip p * in 
( |3a| ) represents the corresponding pump depletion. 

It is often useful to make the simplification of consider- 
ing a situation where the pump, signal and idler energies 
are all close to the corresponding polariton resonance val- 
ues, and the broadenings are small compared to the Rabi 
splitting Q. Then it is a good approximation to replace 
the polariton response by a single Lorentzian function at 
each k, with strength \x\ 2 , where x is the exciton am- 
plitude (Hopficld coefficient) for the mode. The driving 
terms on the right hand side of Eqs.(^) can similarly be 
approximated by (c/ x)f , where c is the photon ampli- 
tude. Then, Eqs.(g) reduce to 

1 C 
-2-(Wp + i'-fp - ujp) ipp + 2~Kip s ipiip p * = — f p (4a) 



II. PARAMETRIC AMPLIFIER 

In the parametric amplifier, both f p and f s are non- 
zero. It is also helpful to assume that f s -C f p , so that 
ip s and ipi are small, and the pump depletion term in 
Eqs.( fja| , ff"a| ) can be neglected. Consider first the situation 
when the probe, idler and pump satisfy the triple reso- 
nance condition, so Eqs.(^) can be used with ui p = uip, 
ui s = ui® and uii = cj?. Without the pump depletion term, 
these equations are solved by eliminating ip p and ipi using 
(Q) and then ([h]), to get 



x, 



x,, 



-V., 



2 -(u)° + ij s - ui s ) ip s + Kip p 2 ipi* = — fs (4b) 

X s 



2„/. * 



X; 



:(ui° + IJi - UJi) Ipi + Klpp i> 



(4c) 



where uip, ui® and ui® are the polariton resonance frequen- 
cies, and 7 P , 7 S , 7, the corresponding widths. In writing 
the equations in this form, the exciton renormalisation is 



-i c s x s *f s /j s 



l-K 2 ^ Xi 



A'., 



H7s 



(6) 



\f P 



Dividing by the value of ip s without the pump, ie with 
fp = 0, gives the internal gain for the probe: 



(7) 



l-Pp/ll 



where 



la = 



lp \flsli 



K C v 



x„ 



Xs Xi 
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The gain increases from unity at I p = to become sin- 
gular at I p = Iq. This suggests that I represents the 
threshold pump intensity for oscillation, which will in- 
deed be shown to be the case in the next section. 

The previous discussion was limited to the case where 
all the fields are on resonance. If this condition is not 
satisfied, it is still possible to obtain an analytic solution 
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FIG. 1. Dependence of the signal response on the pump 
angle, with a probe at normal incidence. The pump amplitude 
f p — 0.292, which corresponds to I p = 0.757o for this geom- 
etry. Structure parameters are: Rabi splitting Q = 5.0 meV, 
zero detuning, exciton width 7^= 0.25 meV, cavity width 7 C = 
0.25 meV, nonlinearity n = 1. The inset shows the polariton 
dispersion with the pair scattering from pump to signal and 
idler modes. 



when the exciton renormalisation is neglected. Solving 
Eqs.(|) with the same weak probe approximation gives 



-(fi/2)/ s /A s 



(uj c (k s )-uj s ) (LUcjhy -ujj) (n/2)* 



1 -K- 



\U\ 



(9) 



A; 



where A p = (lo x — aj p )(uj c (k p ) — ui p ) — (f2/2) 2 , and A 
are similarly defined. 

Fig|l] shows the signal response at normal incidence, 
calculated using Eq.(0), for different pump angles. The 
pump energy is varied to be on the polariton resonance 
for each angle. The spectra show the two polariton fea- 
tures at normal incidence, with clear gain on the lower 
branch for pump angles in the region of 16°: for this 
pump angle, the pair scattering resonance condition is 
satisfied when the probe angle is 0°. 

In Figj2[ the signal response is mapped out as a func- 
tion of energy and angle with the pump kept on reso- 
nance at 16°. There is an enhanced probe response when 
cither the signal or the idler is on resonance - these sin- 
gle resonance energies are indicated by the dashed lines 
on the figure. The response is much stronger for the sig- 
nal resonance, because the probe couples to the signal 
directly, but to the idler only via parametric scattering. 
The strongest response occurs at the double resonances, 
where the dashed lines intersect (at 0°, 16° and 33.5°), 
but there is a long segment of the dispersion very close 
to double resonance, where the gain remains high. 

III. PARAMETRIC OSCILLATOR 

In the parametric oscillator regime, there is no probe, 
so f s = 0, but solutions can still be found with finite 



FIG. 2. Contour plot of the signal response \ip s \ 2 as a func- 
tion of probe energy and incidence angle. The contours are 
logarithmically spaced, and the light regions of background 
shading indicate the highest intensities. Dashed lines show 
the single resonance conditions for signal and idler. The pum 
is on resonance at 16°, and other parameters are as in 
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signal and idler fields. Again, it is simplest to start with 
the triply resonant case. Focusing on Eqs. (|4b| , ^ ), taking 
the complex conjugate of one of them, and treating ip p as 
a parameter, there are non-zero solutions for ip s , ipi only 
if the determinant of the coefficients is zero, that is: 



K \lp p 



X, Xi 



(10) 



The physical interpretation of this condition is obvious - 
for a steady state solution, the generation rate of polari- 
tons in the signal and idler directions, on the left hand 
side, must equal the (geometric) mean of the loss rates, 
y'jsji ■ This is only possible with the value of ip p in 
Eq.(|l0|). For a given external driving field f p , the re- 
quired value of ip p is attained by the depletion of the 
pump polariton field due to the stimulated scattering 
term. 

The resulting signal intensity is calculated by using 



Eq.(^b|) to write ipi in terms of ip. 
( ffa| ) to obtain 



then substituting in 



\ip s 



\f P \ 



(ii) 



where \ip p \ is now just a constant given by Eq.(|10|). Since 
the emitted signal intensity, I Sl is proportional to IV's) 2 , 
this relationship is of the form 



x 



(12) 



where Iq is the same threshold intensity as in Eq.(|8|). Of 
course I s > 0, so this solution only exists when I p > Jo- 
The treatment can be extended to the case where the 
signal direction is such that the signal and idler are not 
both on resonance, that is the mismatch A = 2lu p — uj® — 
uji 7^ 0. The steady state condition, Eq. (|l0|) becomes 



k|</v|VsIN = [{u° s + ij s - uj s )(uj9 - iji - , (13) 

with, once again, loi = 2lo p — uj s . A solution is only 
possible if the right hand side of Eq.(p"3|) is real, which 
requires U) a — = Aj s / (j s + 7,: ) , and correspondingly 
Wi — ujf = Aji/^g + 74). The physical significance of 
this requirement can be seen by looking at ip s and ipf. 
for the allowed value of lu s , -f s \^ s \ 2 /\x s \ 2 = li\ipi\ 2 /\Xi\ , 
that is, the loss rates through the signal and idler modes 
are identical. This is the Manley-Rowe relation for the 
parametric oscillator. 

Continuing the solution for the signal intensity, for res- 
onant pumping the effect of the finite mismatch A is to 
shift the threshold, so 



7 (A) = / (0) 1 + 



(7s + 7») 5 



(14) 



where -Zo(O) is the value given in Eq.(g). The threshold 
is lowest when A = and pump, signal and idler are all 
on resonance. 

These results show that the signal intensity is deter- 
mined by the pump depletion, which produces the value 
of |?/>p| 2 required by Eqs.(|l||||). For I p > I , \ip p \ 2 re- 
mains unchanged, just as the population inversion in a 
conventional laser is clamped at its threshold level. Since 
the actual value of IV^I 2 is unique to a particular pair of 
signal and idler directions, in equilibrium there can only 
be one finite signal amplitude. It is easy to see what will 
happen in an out of equilibrium situation when there is 
more than one signal. A signal whose loss rate exceeds 
its generation rate will decay, while one for which the 
generation rate exceeds the loss will grow. So, in a pro- 
cess akin to mode selection in a laser, the signal with the 
lowest loss rate will dominate, depleting the pump until 
only it survives. 

IV. DISCUSSION 

This model of the parametric oscillator makes two sim- 
ple predictions which can be checked against experiment: 
the \TTp power dependence in Eq.(n2h, and the clamp- 
ing of the pump polariton amplitude, ip p , to the value 
given in Eq.(|H^). The latter effect should be observable 
as a saturation, above threshold, of the blue shift of the 
polariton dispersion. Below the threshold the shift is 
roughly linear with pump power, and it saturates at the 
value given in Eq.(|) with I p = Iq. For the signal mode 
the saturated shift is thus Soj® « 2^/7^7^ |x s |/|Xj|. Of 
course, there are other effects, not included in the model, 
which can cause energy shifts. These include the exci- 
ton renormalisation due to the signal field, and at higher 
powers the break down of the strong coupling regime. 
The former effect can be estimated and is small: the sig- 
nal field will give an energy shift of k|x s | 2 |t/; s | 2 , which is 
about 10% of the saturated Sui® when L, = 2J , using the 
same structural parameters as in FignT. 
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FIG. 3. Theoretical fits to the experimental data (points) 
for (a) signal power and (b) blue shift from Ref. || 



Fig.|| shows how these predictions compare with the 
experimental results in Ref. H. The signal power is fit- 
ted to the form I s cx yJTp~ — ^JIq of Eq. ( |i"2| ) . Although a 
good fit is obtained, it is difficult to distinguish the ^/Jp~ 
behaviour from a simple linear increase with the limited 
range of data above threshold. The support for the model 
provided by the blue shift data is rather better: the ex- 
periments show a clear saturation of the blue shift, and 
the saturation shift of just under 1 meV agrees very well 
with the prediction of 8lo° s k 0.97meV obtained using the 
experimental widths 7 S = 0.57meV and 7, = 0.80mcV. 

The predicted form of the gain in Eq.(^) is not found 
in the pump-probe experiments of Ref. |l|, where the gain 
increases exponentially with I p . However, the gain mech- 
anism in the ultrafast measurements differs from the 
CW parametric amplification discussed here. The peak 
pump powers in the experiments are orders of magnitude 
greater than the threshold Iq, but oscillator behaviour is 
not seen because there is insufficient time for the signal 
to build up just from the incoherent photoluminescence. 
However, the probe pulse seeds the signal, starting an 
an exponential growth which lasts as long as the pump. 
This gives the observed exponential dependence on I p . 
Although the CW gain below threshold predicted here 
has not been seen experimentally, it should be observable 
in structures which show parametric oscillator behaviour. 

To conclude, it has been shown that the parametric 
oscillator behaviour recently observed in microcavitics 
can be well explained by a simple classical treatment of 
a non-linear exciton strongly coupled to a cavity mode. 
The main non-classical effects that are missing from this 
model are the incoherent photolumincsccnca, discussed 
theoretically in a recent paper by Ciuti et al.p and quan- 
tum statistical effects in the oscillator behaviour around 
the threshold. It is also possible that at higher excitation 
powers the signal intensity will become large enough for 
the appearance of new physics due to polariton-polariton 
interactions in the coherent state. 
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